The optimal inhomogeneity for superconductivity 
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We study the effect of nonuniform transverse couplings on a quasi-one dimensional superconduc- 
tor. We show that inhomogeneous couplings quite generally increase the superconducting (pairing) 
gap relative to the uniform system, but that beyond an "optimal" degree of inhomogeneity, they lead 
to a suppression of the tendency to phase coherence. The optimal conditions for superconductivity 
are derived. We also show that a delocalized, spin-gapped phase is stable against weak disorder in 
a four-leg-ladder with moderate repulsive interactions. 
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A number of experiments [1] including neutron scatter- 
ing, angular resolved photoemission, and scanning tun- 
nelling microscopy [2] suggest that many high-T c su- 
perconductors show a highly inhomogeneous electronic 
structure which often has a quasi-one dimensional nature 
- "stripes." It is still unclear whether inhomogeneities, 
and stripes in particular, are an essential feature of high- 
Tc superconductivity. 

There is compelling theoretical evidence that quasi-one 
dimensional systems, such as two-leg Hubbard or t — J 
ladders (2LLs), have a strong tendency towards the for- 
mation of a spin gap and substantial superconducting 
(SC) pair-field correlations [3]. However, the gap size 
and the tendency to superconductivity tends to decrease 
rapidly with width for ladders with more legs. [4, 5] Due 
to the phenomenological similarity between 2LLs and 
two-dimensional cuprates, it has been suggested that the 
pairing mechanism for superconductivity in cuprates may 
have a quasi-one-dimensional origin [6]. However, strong 
pairing is often acompanied by a small superfluid stiff- 
ness, especially in quasi-one-dimensional systems, so a 
large spin gap does not necessarily imply a high-T c [6- 
8]. In other words, while local pairing develops easily 
in a 2LL, it is hard for the phase of the pairs on dif- 
ferent weakly coupled 2LLs to develop coherence. An- 
other problem with quasi-one dimensional superconduc- 
tors is their extreme sensitivity to disorder [9]. Since 
weak disorder does not destroy superconductivity in two 
dimensions, it is intuitively clear that this sensitivity is 
mitigated if one increases the number of coupled chains. 
Taken together, these observations suggest there exists an 
intermediate "optimal" degree of inhomogeneity [4, 10] 
which maximizes T c . 

In this Letter, we analyze the relation between inhomo- 
geneity and superconductivity in "microscopic" inhomo- 
geneous multi-leg ladder Hubbard models in which the 
SC gap arises solely as a result of repulsive interactions 
between electrons. This allows us to make quantitative 
estimates (at least for weak coupling) of the optimal de- 
gree of inhomogeneity without ad hoc assumptions. In 
addition, we show that for a 4- leg ladder (4LL), there is 
a broad range of repulsive interactions for which disorder 



is irrelevant in the renormalization-group sense. 

Specifically, we have carried out a weak-coupling 
renormalization-group (RG) analysis of a model (Eq. 1) 
of two 2LLs with a repulsive on-site interaction U cou- 
pled via an interladder hopping t and with an on-site en- 
ergy offset e between the 2LL. We consider t in the range 
1 > t > 0; in the units we have adopted t = 1 and e = 
corresponds to a homogeneous 4LL, while at t = there 
are two (inequivalent) decoupled 2LLs. For simplicity, 
in order to have a unique parameter which determines 
the degree of inhomogeneity, e is taken to depend on t 
as e = eo(l — i 2 ). Moreover, we have considered both a 
flat ladder, with open boundary conditions in the direc- 
tion transverse to the chains (OBC), and of a cylindrical 
ladder, with periodic boundary conditions (PBC). 

Our results, valid qualitatively in a range electron con- 
centrations per site v near but not equal to 1, are pre- 
sented in Fig. 1, and can be summarized as follows: In 
the parameter range that we have considered, there are 4 
bands which cross the Fermi energy, so there are poten- 
tially 4 distinct gapless charge (p) and spin (s) modes. 
Upon bosonization, these modes are represented by the 
collective bosonic fields, <j> pM and <j) Sta where a = 1 — 4 is a 
band index [11]. Phases (that is to say fixed points of the 
RG flows) are labelled CnSm according [5] to the number 
n and m of charge and spin modes that remain gapless 
in the presence of interactions. For the entire range of 
t, we will show that the 4LL is in the maximally gapped 
CI SO phase in which all 4 spin modes and the 3 "relative" 
charge modes are gapped [12]. (So long as v is irrational, 
the total charge mode (f> p+ = (1/2) <j> Pta , is always 
gapless[ll].) This phase is characterized by slowly de- 
caying power law SC correlations, and can be considered 
as the "SC" phase of the 4LL. For a range of t, a sin- 
gle, strong-coupling fixed point governs the physics, and 
hence all the gaps are "comparable" [13], in the sense that 
they all have the same singular dependence on the inter- 
action strength - this single scale behavior is remeniscent 
of an anisotropic BCS superconductor. In other ranges, 
the RG flows pass close to an initial strong coupling fixed 
point before finally reaching the C1S0 fixed point. Here, 
there are two distinct gap scales, a "pseudo-gap" scale, 
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Ape with magnitude governed by the flow to the ini- 
tial strong-coupling fixed point, and a lower scale, which 
we call the superconducting gap Asc, governed by the 
second segment of the flow. 




FIG. 1: Logarithm of Asc and A PG for a 4LL with OBC 
(solid line) and PBC (dashed) for eo = 1, v = 1.2, and U = 1. 
We also indicate the positions of t op t, t* , and t** for OBC 
(upper row) and PBC (lower). 

The crucial result (see Fig. 1) is that Asc is not a 
monotonic function of t, and, in particular, it has an ab- 
solute maximum at a certain t — t opt < 1, i. e. for 
an inhomogeneous 4LL. We find three regimes of t: For 
t** > t > t* , all the gaps are comparable. In contrast, 
two distinct gap scales are found for 1 > t > t** and 
t* > t > 0. The two-gap behavior at small t can be un- 
derstood as follows: For t = 0, one of the 2LLs has an 
exponentially larger superconducting gap than the other 
due to the different Fermi velocities. For small but non- 
zero t, Ape is continuously connected to this gap, i.e. 
it is the gap associated with the pair of bands predom- 
inantly associated with the more SC 2LL. However, for 
finite t, smaller gaps of magnitude Asc are induced in 
the remaining bands by a generalized version of the prox- 
imity effect [6]. 

a. Inhomogeneous 4-LL in weak coupling: We con- 
sider a Hubbard Hamiltonian for two coupled 2LLs 
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H = -H E E H,y,*c*+i,y,° + H.C.] (1) 
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~ E tj ~ [ C lx<J C x,2.a + 4,3,<t C x,4,<7 + H.C.] 
~ E [ tC l,2,<T C 3,x,<? + <'4,4,<r c «>l^ + H ' C -] 
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where, c\ creates an electron with spin a on chain 



at position x along the chain and h x , y — X) 



'/ - . jo C x,y,<T C x,y," 



We choose units in which the hopping matrix elements 
along the ladders, ty = 2. We set t± — 1 along the rungs 
of each 2LL, and a different coupling t < t± along the 



rungs connecting the two 2LLs where t' = corresponds 
to OBC and t' = t to PBC. In addition, one of the two 
ladders is shifted in energy by an amount e = eo(l — 
t 2 ). The electrons interact via a weak, on-site Hubbard 
interaction U -C 27rtii. We fix the number of electrons 
per site, to be close to, but not equal to one. 

b. First stage renormalization: It is by now a 
straightforward procedure [5] to derive the weak-coupling 
RG equations for the large number of distinct running 
couplings gi that define the various low-energy two- 
particle scattering processes. These equations always 
have the form 



dg./dT^A'g^n + Oig 3 ), 



(2) 



where r is the log of the bandwidth (cutoff), and Aj t 
are coefficients which depend on the details of the non- 
interacting system under consideration, i.e. in the 
present case they depend on eoj t, and v. The other 
inputs to the calculation are the initial values gi(r = 
0) = XjU <C 1, where the dimensionless constants Xj also 
depend on eoj t, and v. Information about the charac- 
teristic emergent energy scales is obtained by integrating 
these equations (typically numerically) until one or more 
coupling grows to be of order 1. 

To be precise, a characteristic gap scale, 



A 



exp[ 



exp 



U 



O (log CO 



(3) 



is obtained from the value of r = t* = a/U at which 
some set of couplings diverge. (The 0(\ogU) correc- 
tion represents the effect of the neglected cubic terms 
in powers of g.) For r near r*, there are generally a 
set of strongly divergent couplings, gi ~ Gj(r* — t) _1 , 
which all grow to be of order 1 when (r* — r) ~ 1, 
and so define the new "strong-coupling fixed point." 
There is also a set of weakly or non-divergent couplings, 
gi ~ GiU Ji (T* — t)^ 1 k > with 7$ > 0, which remain 
small, gt - G t U^, when (r* - r) ~ 1 [5, 14]. What 
modes are gapped at this strong-coupling fixed point is 
then determined by bosonizing the model. The order 1 
couplings are considered to be large and thus pin (gap) 
the appropriate modes. As we will discuss below, it may 
still be necessary to do further analysis to address the fate 
of the remaining gapless modes at this strong-coupling 
fixed point, and to account for the effects of the weak 
residual interactions between them [6]b. The coefficient 
a in Eq. (3) , which is the principle result of this initial 
analysis, clearly depends on the the values of t eo, and v. 

Since this analysis straightforwardly extends the pre- 
vious treatments of the 2LL and 4LL [5], we will not, 
here, present the details of the calculation. Some rep- 
resentative results, for parameters listed in the captions, 
are shown in Fig. 1. For homogeneous 2iVLLs, a in- 
creases with the number of chains 2N, as pointed out in 
Ref. [5], so that the gaps decrease exponentially with N. 
In the present case, a(t) must interpolate between the 
relatively large value for a homogeneous 4LL when t = 1 
and the smaller 2LL value when t = (See Table I.) . 
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By integrating the one-loop RG equations, we found, to 
our surprise, that even the dependence of a on t is not 
monotonic; a(t) first increases with decreasing t until it 
reaches a maximum value at t = t min , and then decreases 
(i. e. the gaps increase exponentially) as t decreases fur- 
ther. 

For t* < t < t** (See Table I.) all the couplings re- 
sponsible for gapping the 3 charge and 4 spin modes di- 
verge in proportion to (r* — t) _1 at r = t* , so that 
all gaps at this C1S0 fixed point are "comparable" [13]. 
On the other hand, for t* > t > or 1 > t > t**, 
only the couplings responsible for the gaps in the spin 
and relative charge modes within two bands diverge like 
(t* — t) _1 . More specifically, at the resulting (C3S2) 
fixed point, Ape = A characterizes the pinning of the 
spin fields [11] <f> Si i, <p s ^ associated with two of the bands 
(labelled, for simplicity, 1 and 2), and of the relative "su- 
perconducting phase," 0^,(1-2) = (9 p a — 0p t 2)/V%, where 
Qp.a is the field dual to (f> p . a - Behavior of this sort has 
sometimes been interpreted [5, 14] as indicating a tran- 
sition to a phase with additional gapless modes, i. e., 
in this case, from a C1S0 phase for t** > t > t* to a 
C3S2 phase for t < t* and t > t** . However, we shall 
see[6]b that there is a crossover, but no phase transition, 
at t = t*. 

c. Second-stage renormalization: For t < t* , the 
C3S2 strong-coupling fixed point to which the weak- 
coupling flows have carried us can be thought of as de- 
scribing a ID superconductor in bands Oa = 1,2 and 
two ungapped Luttinger liquids (LL) corresponding to 
bands bp = 3,4. The fixed-point Hamiltonian thus con- 
sists of 3 gapless charge and 2 gapless spin modes with a 
renormalized ultra-violet cutoff Ape- However, several 
residual interactions are left at the end of the first stage 
of renormalization, and it is necessary to carry out a per- 
turbative stability analysis of the C3S2 fixed point with 
respect to these interactions [6]. In the weak-coupling 
limit, it turns out that the four singlet Josephson cou- 
plings (Jb A ,b B ) between the bands bA and bp, i.e. the 
amplitudes to scatter a zero-momentum pair from bA to 
bp, are the most relevant perturbations and make the 
C3S2 fixed point unstable. Other interactions, such as 
the interband or intraband scattering processes in the 
bands B, are either irrelevant or marginal in the 17 — *■ 
limit. 

This can be seen by evaluating the scaling dimension of 
the Jb A ,b B {dj, approximately equal for both bp), which 
can be estimated again using bosonization. Specifically, 
one replaces the gapped fields of the two bands 6 a with 
their expectation values, while all other fields can be 
taken as nearly free, i. e. their (spin and charge) LL 
exponents K can be approximately set to 1, since the 
remaining couplings are small. This procedure (which is 
standard [6]) leads to the estimate 



dj w 1 + 1/(4 K pA ) 



(4) 



where K p a is the charge Luttinger exponent of the to- 
tal charge-mode of bands bA- Since [15] dj w 5/4 < 2, 



Jb A -,b B are relevant perturbations, and diverge at lower 
energy scales within a strong-coupling RG expansion 
about the C3S2 fixed point, as anticipated. Again, 
one can identify the energy scale at which one or more 
Jb A ,b B become of order 1 with the freezing of the as- 
sociated modes and with the opening of corresponding 
gaps. Specifically, Jb A ,b B open a gap in the spin modes 
(<t> s ,b B get locked) as well as in the relative charge modes 
(locking the relative SC phases (0 p a +0 P ,2 — 26 p ,b B )/V6) 
resulting in a C1S0 phase. Since all Jb A .b B are compara- 
ble [13] and have the same scaling dimension [Eq. (4) ], 
all these gaps are comparable as well. These smaller gaps 
are thus identified as Asc of the 4LL. In this C1S0 phase, 
only the global charge mode <j) p+ remains ungapped. The 
criterion discussed above yields 



^sc 



fl FG Jab 



(5) 



Here, J\ B is the value of the residual Josephson cou- 
plings J~b A ,b B left at the end of the first stage of renor- 
malization, i. c. when (t* — t) ~ 1. The numerical 
asymptotic analysis around the singular point r* yields 
a weak divergence of the J~b A ,b B with a unique exponent 
Iab ~ 0.06. Thus, the 3~b A ,b B attain a value 



Jb 



b A ,b E 



U 1AB G AB 



(6) 



which is small but finite for small U. (It is interesting 
that jab ~ 0.06 is independent of t (See Table I.) , 
although clearly Gab vanishes for t — > 0.) 
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tmin 
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a{t = 1) 


JAB 


r-(2) 

AB 


OBC 


0.65 


0.91 


0.81 


20.1 


16.8 


35.1 


0.06 


0.05 


PBC 


0.92 


1.0 


0.91 


20.1 


38.9 


41.8 


0.06 


0.01 



TABLE I: Values of the parameters discussed in the text for 
a 4LL with OBC and PBC, and for v = 1.2, and e = 1. 

The optimal t is not given by t = t* . Although the ratio 
A sc /Ap G is a rapidly decreasing function of decreasing 
t for t < t* , Asc itself still, in general, initially increases 
due to the exponential increase of Ape- Anticipating the 
fact that the optimum value for t 2 is of the order U, we 
can determine it by expanding a for small t as a(t) = 
(Xo + cvit 2 + 0(t 4 ), with ai given in Table I. Similarly, 

Gab = G% t 2 + C(t 4 ). Then from Eq. (5) and from 
the expression for J\ B 



A sc (i)~ A PG (t = 0) 



JJ1A 



u AB t 



- ai t 2 /U 



(7) 



where xj = 1/(2 — dj). It is now straightforward to 
determine the value of t which maximizes Asc- 



t z - 



U 



{2-dj)a_ 



+ 0{U 2 ) 



(8) 



d. Quantitative considerations: The various quanti- 
ties describing the t dependence of the gaps discussed 
here are reported in Table I for v = 1.2 and eo = 1. The 
curves in Fig. 1 have been computed as follows. We have 
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taken initial values of the coupling constants correspond- 
ing to U = 1 and have integrated Eq. (2) numerically for 
different values of t until the largest of the couplings (in 
modulus) has reached the value 1. This determines r* (t), 
and where needed, J\ B . The gaps are then evaluated by 
using Eq. (3) , and Eq. (5) with dj = 5/4, respectively. 

A different situation occurs when eo = 0. In this case, 
the wave functions of the U = t = part of the Hamil- 
tonian are equally distributed on the two 2LL for any 
nonvanishing t. As a consequence, the dependence of 
the gaps on t obtained by the method described above 
is very weak (for PBC there is no dependence at all). In 
this case, one should use a different method in which cou- 
plings with crystal momentum not exactly conserved at 
the Fermi points have to be included in the first stage 
of the renormalization. Qualitatively, the results are 
similar to the e ^ case, i. e., the SC gap increases 
with increasing inhomogeneity, there is an optimal value 
t = t op t, an d a "pseudogap" behavior. However, in this 
case t op t is much smaller and of the order of the 2LL gap 
A PG (t = 0) [16]. 

e. The effects of disorder. Orignac and Gia- 
marchi [9] have shown that disorder is a relevant per- 
turbation (i.e. leads to localization) in the C1S0 phase 
of the 2LL unless there are strong enough attractive in- 
teractions such that the Luttinger exponent, K p+ , asso- 
ciated with the gapless total charge mode is greater than 
3/2 [15]. We now consider the effect of a single-particle 
disorder potential in the C1S0 phase of the 4LL. 

Neglecting forward scattering terms, which turn out 
not too be crucial here, one can write this potential as 



r 4 

J U Uf 



Vb,b'(x)e 



H.C. (9) 



Here, the operator Vb,b' ( x ) scatters an electron from band 
b to band b' and from left to right, v^y is straightfor- 
wardly expressed in terms of the bosonic fields. Since in 
the C1S0 phase, all spin-fluctuations are gapped, to study 
the low energy consequences of weak disorder we can re- 
place all factors of the form e la ^ s - b by its constant expec- 
tation value. The remaining dependence on the charge 
and dual spin fields is 



Vb,b>{x)~— e ^ cos 

lira 



7 s,b 



7s,b' 



V2 



+ V 



(10) 

where the x dependence of the bosonic fields is left im- 
plicit, and r) is a x-independent phase. The £b,b> {%) 
are the random potentials which, for the sake of defi- 
nitencss, can be taken to obey a Gaussian distribution 
£.b u b 2 (x)£b 3 ,b 4 (y)* = W blM \ b3M D5{x - y), where D is 
the disorder strength (proportional to the inverse of the 
mean- free path) , and W are numbers of order unity that 
depend upon details of the band structure. The precise 
form of the W as well as the value of T] are not important 
for the following analysis. 

Now we get to the nub of the problem. The only low 
energy fluctuations in the CI SO phase involve the total 



charge mode, 4> p +- Since the spin- fields are pinned, the 
dual spin fields are wildly fluctuating, so that any opera- 
tor that depends on e la0s < b has exponentially falling corre- 
lations, and is hence irrelevant in the RG sense. Similarly, 
since the relative superconducting phases, 9 Pi b — 9 P ,b>, of 
the various bands are pinned, any operator with a fac- 
tor e M ^ p ' 6_l *p. 6 '] is similarly irrelevant. Therefore, both 
terms in Eq. (10) are irrelevant at low temperatures. 

This does not mean disorder is irrelevant. As discussed 
in Refs. 17, 18, additional terms get generated in early 
stages of the RG flow as the gapped fields get integrated 
out. It is straightforward to show that the most relevant 
term is obtained only at fourth order in v, and has the 
form 



H, 



eff 



dx Zeff(x) eiV S. 

2ira 



H.C. 



(11) 



where £ e ff is an effective disorder potential, which 
is again Gaussian distributed: £,eff(x)£ e ff(y)* = 
D eff S(x - y), where D eff ~ D(D/v 2 F ) 3 . 

For weak enough bare disorder [19], and for energies 
below the SC gap Asc, the RG flow for the renormalizcd 
disorder strength D e ff associated with Eq. (11) can be 
easily derived in the usual way [9]. The RG equation 
reads 



d D 



eff 



d 



= D eff (3-AK p+ ) 



(12) 



i. e. disorder is irrelevant for K p+ > 3/4. This results is 
very important, as it signals the presence of a localized- 
delocalized transition in a quasi-one-dimensional system 
for purely repulsive interactions [15]. Moreover, the fact 
that delocalization is present in weak coupling, i. c. 
within the range of validity of the present RG proce- 
dure, makes the 4LL one of few low-dimensional mod- 
els in which one can show the occurrence of delocal- 
ized states in a controlled manner. This result is valid 
provided the disorder does not destroy the gaps, i. c. 
D < v F A sc [19]. 

/. Comments and Speculations. The present results 
were derived for small U, where all the emergent energy 
scales are exponentially small. In many cases, the same 
physics, and in particular the existence of an optimal de- 
gree of inhomogeneity, can be shown[6, 16] to apply in 
the strong-coupling limit, as well, where the gap scales 
are a significant fraction of the exchange interaction, J. 
However, quantitative results in this regime generally re- 
quire some input from numerical experiments. 

We note that the present results for the 4LL are qual- 
itatively similar to what happens in high-T c supercon- 
ductors such as La2- K Sr :E Cu04 [20, 21]. In the over- 
doped region the system is relatively homogeneous, both 
the pairing and SC scales are relatively small and roughly 
equal. With decreasing doped hole concentratin, x, stripe 
correlations become more prominent and the interstripe 
distance increases, i. e. the system becomes more effec- 
tively inhomogencous. At the same time, the pseudo-gap 
scale increases monotonically with decreasing x, while 
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the SC T c has a maximum at a "optimum" x and then 
decreases. Besides High-T c superconductors, the effects 
studied here coutd be relevant for quasi-one-dimensional 
systems, such as quantum wires or carbon nanotubes, in 
which more than one channel cross the Fermi surface. 
Finally, this idea could, possibly, give direction to the 
search for better superconductors [8]. 
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